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Abstract 

We show that in order to calculate correctly the spin current carried by a quasiparticle in an 
electron liquid one must use an effective "spin mass" m s , that is larger than both the band mass, 
rrib, which determines the charge current, and the quasiparticle effective mass m*, which determines 
the heat capacity. We present microscopic calculations of m s in a paramagnetic electron liquid in 
three and two dimensions, showing that the mass enhancement m s /mb can be a very significant 
effect. 
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The recent explosion of interest in metal and semiconductor spintronics |2| has brought 
into sharp focus the basic problem of calculating the spin current carried by a nonequilib- 
rium electronic system. The standard approach is to solve the Boltzmann equation for the 
nonequilibrium distribution function; but this is not sufficient when many-body effects due 
to electron-electron interactions need to be taken into account. In fact, electronic correla- 
tions are particularly strong in low-dimensional systems, such as magnetic semiconductor 
films and wires, which are currently being considered for the realization of spin transis- 
tors [3]. One might hope to take care of the many-body effects by solving, instead of the 
Boltzmann equation, the Landau-Silin transport equation for quasiparticles But even 
this is not sufficient, since the transport equation per se does not tell us how to connect the 
quasiparticle distribution function to the spin-current. The key question, which seems to 
have been overlooked so far in the growing literature on spin transport p, is also a very 
basic one, namely, what is the spin-current carried by a single quasiparticle of momentum 
p and spin a? Without knowing the answer to this question it is not possible to calculate 
the spin current from first-principles. In this paper we show that, in order to calculate the 
spin current correctly, one must recognize that the effective spin mass m s , which determines 
the relation between the spin current and the quasiparticle momentum, is neither the band 
mass mb (which controlls the charge current), nor the quasiparticle mass m* (which controls 
the heat capacity), but rather a new many-body quantity, controlled by spin correlations. 
Our calculations show that the spin mass, in spite of uncertainties due to the approximate 
character of the many-body theory, can be considerably larger than the bare band mass in a 
two-dimensional electron gas (by contrast, the quasiparticle effective mass is typically very 
close to the band mass). Hence, the spin mass will have to be taken into account whenever 
a quantitative comparison between theory and experiment is desired. 

Let us begin by describing the physical origin of the spin mass. The spin current, j s = 
jj — jj, is defined as the difference of the up-spin and down-spin currents, and which 
in turn are defined as the expectation values of the operators 



in the appropriate nonequilibrium state. Here a — 1 for \ spins and a = — 1 for J, spins, 
nib is the bare band mass, pi is the canonical momentum operator of the z-th electron, 
a z ,i is the Pauli matrix of the ^-component of the spin of the z-th electron, 1+(T ^' Z is the 




(1) 
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projector on the cr-spin component of the i-th electron, and N is the number of electrons. 
Let us consider a many-body state, denoted by [per), which contains a single quasiparticle 
of momentum p and spin a. This state carries a total current j = whether or not 
interactions are taken into account. The reason why this is so is simply that the state \p&), 
which contains a quasiparticle of momentum p and spin a, is an eigenstate of the current 
operator j = J^ i=1 ^- with eigenvalue As a consequence, the current density associated 
with the distribution n a (r,p,t) is given by [y] 

J(r, t) = Y] —n a (f,p, t) . (2) 

per 

The difficulty in calculating the spin current arises from the fact that the state \pa) is 
not an eigenstate of j\ or 3 \'- thus, we cannot automatically say that in this state j a = 

and j (j = 0, even though these expectation values would be consistent with the total value 

of j. All we can say, a priori, is that the expectation values of and j± in the state \pa) 
must be proportional to p and add up to Thus, we write 

{w\jr\pcr) = a Ta — , (3) 

where a Ta is a 2 x 2 matrix whose columns add up to 1, so that the total current is 
Notice that in a paramagnetic system = an, and, therefore an = a if. for simplicity's 
sake, we will focus on just this case from now on. The above Eq. (JHJ) implies that the spin 
current carried by an up-spin quasiparticle of momentum p is 

J s (pT) = (<*tt = — > ( 4 ) 

1 1 Of) lll s 

and, similarly, the spin-current carried by a down-spin quasiparticle is 

3s(pi) = (a n -an) — = — — , (5) 
rrib m s 

since att = „ u . The S e e q ua tl0 n S define a sp in m ass m wfiiefi eon trols tfie spin ee rre nt in 
much the same way as m& controls the charge current p|- 

Combining Eqs. (0J) and (0) we see that the correct expression for the spin current density 
carried by a nonequilibrium quasiparticle distribution n a (r,p,t) is 

3s if, t) = — K ( f > P> *) - n l f, P,t)\ ■ (6) 

- Wis 
V 
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FIG. 1: Spin-momentum separation in a Fermi liquid: the momentum of an up-spin quasiparticle 
resides, in part, with down-spin particles. The solid lines represent the Fermi surfaces and the 
dashed lines indicated the collective shifts in the momentum distribution of up- and down- spins. 



It is clear that m s must be larger than mj since and are positive numbers that 
add up to 1, implying that «n — < 1. The positivity of «n and can be intuitively 
grasped by considering the physical picture illustrated in Figure 1. We start from an exact 
eigenstate of the noninteracting system with full Fermi spheres of up- and down-spins and 
an additional single particle of momentum p and spin f out of the Fermi sphere. In this 
state ji = ^- and j± = 0. The quasiparticle state is now obtained by slowly turning on the 
electron-electron interaction. The total momentum and spin do not change in the process, 
but some momentum is transferred from the up- to the down-spin component of the liquid: 
one may say that the up-spin quasiparticle drags along some down-spin electrons as part of 
its "screening cloud". As a result, the magnitude of is smaller than by an amount 
equal to The magnitude of the spin current is a fortiori smaller than which implies 
m s > rrib. Notice that the "spin- momentum separation" described above is entirely due to 
correlations between electrons of opposite spin orientation. Interactions between same-spin 
electrons do not contribute to this effect. 

Having thus clarified the general concept of the spin mass we now proceed to (1) relate 
m s to the quasiparticle effective mass and the Landau Fermi liquid parameters, (2) relate m s 
to the small wave vector and low frequency limit of the spin local field factor G_(q, u), and 
(3) present approximate microscopic calculations of m s in a paramagnetic electron liquid in 
three and two dimensions. 

Let us start from the quasiparticle state \pa) and apply to it the unitary transformation 



U = exp 



1+T<T Z 



, which boosts the momenta of the r-spin electrons by q T . By 



applying U to the fundamental hamiltonian of the electron liquid one can straightforwardly 
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show that the change in energy of any state, to first order in q a , is 



(7) 



On the other hand, for the quasiparticle state under consideration, we know that j 7 
Substituting this into Eq. (JJJ) we get 



AE = Y1 



r 



(8) 



The energy change under this transformation can also be calculated with the help of the 
Landau theory of Fermi liquids. There are two contributions: one from the boost in the 
momentum of the quasiparticle, and the other from the collective displacement of the Fermi 
surfaces by q T . A standard calculation gives 



V 



S„n- — 



(9) 



where n 0tT (p) = OQop — p) is the momentum distribution in the ground state and is the 
Fermi momentum. 

Comparing Eqs. (jS)) and Q, we arrive at the identifications 

F Tt" 



Q!tT 



m* 



1 + 



2d 



m b F{ 
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Qttj = — -^t- , (10) 

u m* 2d K J 

where d is the number of spatial dimensions and F[ T = A r *(0) f ^f^ at f T Pe(cos9) is the 
angular average of the interaction function, weighted with the Legendre polynomial P^(cos#) 
(or just cos £9 in two dimensions) and multiplied by the density of states at the Fermi surface, 
N*(0). Notice that the sum rule ^ T a T(J = 1 is satisfied by virtue of the well known Fermi 
liquid relation ^ 

= I (ii) 

m* 1 + Ff/d ' v ' 

where F^ = | FP + (-)Fj^ are the standard dimensionless Landau parameters defined, 
for example, in Ref. Q. The spin mass, on the other hand, is given by (see Eq. (jl])) 

1 



m s 
m* 



1 + Ff/d ' 



(12) 
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showing that the relation of m s to m* is to the spin-channel what the relation of mj to m* 
is to the density channel. 

It should be noted that the spin current density obtained from Eq. (JHJ) satisfies the 
continuity equation 

^-n s (f,t) + V-j s (r,t) =0 , (13) 

where n s (f,t) = n^(f,t) — ni(r,t), is the spin density. Conversely, Eq. (J12)) could have been 
directly obtained from the requirements of charge and spin conservation. 

The microscopic calculation of Landau parameters is notoriously difficult. Diagrammatic 
calculations of F( and Ff in the three-dimensional electron liquid were done by Yasuhara 
and Ousaka 8j, and the calculated parameters, together with the resulting values of m s /mb 
are listed in the upper half of Table for various values of the Wigner-Seitz radius r s . In 
two dimensions the parameters F* and F" were calculated by a variational Quantum Monte 
Carlo method in Ref. [^|. The parameters and the resulting values of m s /mi, are listed in 
the bottom half of Table UJ Notice that the spin mass enhancement in two dimensions is 
considerably higher than in three dimensions. 

In view of the uncertainty in the calculation of the Landau parameters it seems worthwhile 
to attempt another kind of calculation, which does not rely on diagrammatic expansions. 
We first establish the connection between the spin mass and the dynamical local field factor 
in the spin channel. We recall that the dynamical spin susceptibility of an electron liquid is 
usually represented in the form 

where Xo(q, is the noninteracting spin susceptibility (i.e., the Lindhard function), v q is the 
Fourier transform of the Coulomb interaction (= Aire 2 /q 2 in three dimensions and 27re 2 /g 
in two dimensions) and G^(q,u) is the dynamical local field factor in the spin channel. In 
the limit q — > and small, but finite frequency (uo <C ep/H where ep is the Fermi energy), 
Eq. (fTlj) reduces to 

/ \ g^o nq 2 
Xs{q,uj) -> - 



1 + li m n q *v q G-( q ,u) 



U 2 



(15) 
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d r s 1 2 3 4 5 



3 ^ -0.0543 -0.0647 -0.0713 -0.0773 -0.0829 



^ -0.0645 -0.0825 -0.0915 -0.0956 -0.0965 



1.011 1.019 1.022 1.020 1.015 

m b 



2 ^ _o.071 -0.050 -0.015 0.061 



^ -0.096 -0.120 -0.130 -0.136 



1-028 1.080 1.132 1.228 



TABLE I: Landau parameters from Refs. [8j and |9J and spin mass enhancement ^ = i+^a/^ in 
the d-dimensional electron liquid. 



On the other hand, the small-g/finite-u; limit of Xs{q,^) can also be calculated by solving 



the kinetic equation [4] in the presence of slowly varying external fields V a (q,uj). In this 

region collisions are irrelevant, and one gets the spin response 

2 

Sn s {q,u) ~ V s (q,u) , (16) 

m s u z 

where 5n s (q,uj) = 8n^{q,uj) — 8ni(q,uS) and V s (q,uj) = [V^(q,u) — Vi(q,u)]/2. Therefore, 

Xs[q,uj) -> . (17) 

m s uj z 

Comparing the above equation with Eq. ()15j1 leads to the identification 

^ = l+limlim ng2 ^-^" ) . (18) 

The order of the limits is, of course, essential. When uo tends to zero first, G^(q, u) vanishes 
as q d ~ l for q — > 0, so as to yield a finite enhancement of the uniform static spin susceptibility 



In Eq. (HHJ), however, q tends to zero first, and we see that v q G-(q, uj) must go as ^ in order 
to give a finite value of the spin mass. 

The above analysis, combined with the Kramers- Kronig dispersion relation, leads to the 
following relation between the real and the imaginary part of G-(q,uj): 

lim3?eG_(g,a; = lim-P / du/ ,. , 9 V %. ; , 19 

?->0 ^' g^O 7T Jo Uj'{uj' 2 - UJ 2 ) 

where V denotes the principal-part. In the uj —>■ limit, comparison with Eq. (|18j) yields 

^ = 1 + 2^2 r duJ <l\ZrnG M 
nib ™>b 9-*° tt Jo k> 

The quantity lim^o q 2 v q '3mG_(q, uj) was written as a convolution of the response func- 
tions X(ra'(k,uj) in the mode-decoupling theory [lfl, 11, 12, 3| in 3-d in Ref. |l^, (in which 
it was denoted as A(lu), and note that the prefactor becomes i n d dimensions, with 
L the linear size of the system). The results of our calculations of the spin mass from Eq. 
(|20|). with the response functions Xaa'{k,uj) evaluated in the generalized random phase ap- 
proximationfGRPA), are listed in Table ITT] The static local field factors in GRPA are taken 
from Ref. ^] in 3-d, and from Refs. |l?| and [lsl ] in 2-d, respectively. Although there are 
considerable differences between the numbers obtained in different approximations, we see 
that the values of the spin mass obtained by this method are consistently larger than the 
ones listed in Table |U 

In summary, we have shown that the calculation of the spin-current in an electronic system 
is a delicate task: it is not sufficient to include interactions in the transport equation for the 
quasiparticle distribution function: one must also use the correct spin mass to calculate the 
spin-current from the distribution function. We have found that the difference between the 
spin mass and the bare band mass is much larger in two dimensional systems than in three- 
dimensional ones. Although the spin masses calculated in various schemes in 2-d are quite 
different from each other and might be overestimated in some cases due to the limitations 
of the approximations employed, there is no doubt that they all indicate a significant many- 
body effect which is definitely large enough to be observable in the exciting practice of 2-d 
spintronics. We hope that these results will stimulate more accurate calculations of the spin 
mass by quantum Monte Carlo methods. 

We gratefully acknowledge support by NSF grants DMR-0074959 and DMR-0313681 and 
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d r s 1 2 3 4 5 



6 m b 



1.02 1.06 1.11 1.17 1.23 



Ills 

m b 



1.01 1.03 1.03 1.04 1.04 



Ills 



1.15 1.46 1.83 2.21 2.59 



1.77 2.78 4.11 5.36 



TABLE II: Spin mass enhancement calculated from Eq. (|2U[). The first, third, and fourth lines are 



calculated with g(u) = 1 



third moment sum rule in 3-d Il3 
from Refs. (third line) and .1.8] (fourth line) in 2-d 



while the second line includes an empirical correction for the 
15j. The local field factors are taken from Ref. [16( in 3-d, and 



by DOE grant DE-FG02-01ER45897. 
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